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Abstract 

We study jamming in granular mixtures from the novel point of view of extended hydrodynamics. 
Using a hard sphere binary mixture model we predict that a few large grains are expected to get 
caged more effectively in a matrix of small grains compared to a few small grains in a matrix of 
larger ones. A similar effect has been experimentally seen in the context of colloidal mixtures. 
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I. INTRODUCTION 

In recent years granular matter has emerged as an active area of research [l,]. Interest in 
this field has grown as a result of observations coming from a large number of interesting and 
relatively low-tech experiments [2J. In particular, a vibrated granular system^ consisting of 
a large number of macroscopic grains in motion, provides us with a "scaled up" fluid system 
where we can explore the similarities and differences between such a large scale system and 
a microscopic atomic fluid. 

Recent studies in dense granular systems indicate that there is a striking similar- 

ity between the dynamics of granular materials and glassy dynamics in atomic fluids and 
colloids 6|. In this Letter, we exploit this analogy to understand some facets of jamming 
in granular matter. In particular, we study the role of compositional disorder in form- 
ing jammed configurations which slow down the dynamics and eventually result in a state 
of structural arrest characteristic of a glass 0, ls[ 3, [hJ. The observations made here are 
analogous to jamming effects studied in atomic liquid mixtures and colloids 

an 112). 



The behavior of a fluid at large length and time scales is well described by hydrodynamics. 
The hydrodynamic description has been extended to molecular length scales in an extended 
or generalized hydrodynamic description [?], [l3]. In a generalized hydrodynamic framework 
the basic structure of the hydrodynamic equations is retained and the static susceptibili- 
ties and transport coefficients are wave- vector dependent to account for the nontrivial static 
correlations that come into play on molecular scales in a dense liquid. The generalized hydro- 
dynamic description has proved to be very successful in describing the dramatic narrowing 
of the central diffusive peak in a neutron scattering spectrum S(k,u) of a dense liquid [141] . 
Such a slow decay of density fluctuations is a precursor to glassy dynamics in dense liquids. 
A binary liquid mixture turns out to be more effective in forming glasses compared to a 
one component liquid since compositional disorder leads to jammed configurations which 
prevent the system from reaching its global equilibrium crystalline configuration 15]. Gen- 
eralized hydrodynamic studies of dense binary hard sphere mixtures in the context of glass 
transition suggest that at intermediate wavevectors (i.e. equivalently on length scales of the 
order of the average of the molecular diameters of the two species) the density fluctuations of 
the two species emerge as the slowest decaying fluctuations and therefore dominate the slow 
dynamics of the system [?], Jjj]- In a self-consistent mode-coupling theory (MCT), nonlinear 
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couplings of these slowly decaying modes of density fluctuations lead to a glassy state where 
structural relaxation is frozen j^, 9]. 

In this Letter we use, for the first time, generalized hydrodynamic techniques developed in 
the realm of atomic liquids to understand the behavior of jammed configurations in granular 
mixtures. Such an approach to granular matter enables us to take into consideration non- 
trivial static correlations stemming from the granularity or finite sizes of the particles, which 
in turn, influences the high density jamming dynamics of the system. We thus gain new 
insight into granular mixture dynamics and make predictions for future experiments in such 
systems. 

II. EXTENDED HYDRODYNAMIC APPROACH TO A BINARY MIXTURE OF 
GRAINS 

Consider a binary granular mixture of hard spheres of diameters 0\ and cr 2 (<7 2 >(Xi), 
masses ni\ and m 2 , number densities n\ and n 2 and of total packing fraction rj = 
^[n\af + n 2 o"|]. An analysis of the extended hydrodynamic equations of such a binary 
hard sphere mixture suggests that on length scales of the order of the size of the hard sphere 
particles |7| momentum and temperature fluctuations decay very fast and the slow dynamics 
of the system can be well described in terms of the modes of density fluctuations of the two 
species This enables us to confine to a two mode description of the system on length 

scales of the order of the diameters of the hard sphere granular particles at high densities. 
It is convenient to describe the system in terms of the following linear combinations of the 
density fluctuations of the two species- the total mass density fluctuation at wavevector k 

Pk = Pik + P2k 
and the concentration fluctuation at wavevector k 

_ P2 _ fh 

°k - p 2 Plk p 2 P^k 

where p x % and p 2 £ are the mass density fluctuations of species 1 and 2, pi = mini and 
p 2 = nfi2U2 are the equilibrium mass densities of species 1 and 2 and p = pi + p 2 is the total 
equilibrium mass density. The set of Laplace transformed coupled extended hydrodynamic 
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equations of this system is given by: 
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where Xt(^) is the generalized isothermal compressibility, defined in terms of the partial 
static structure factors Sij(k) with i = 1, 2 and j = 1, 2: 

S n (k)S 22 (k) - Sf 2 {k) 
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x 2 S n (k) + xiS 2 2(k) - 2y/xix 2 Si2(k) 
Here Xt ^ s the compressibility of the granular gas in the dilute limit. 7l(/c) is the generalized 
longitudinal viscosity [?}]. 
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Sn(k)S 22 (k) - S 2 2 (k) 

are combinations of partial static structure factors and D(k) is the coefficient of mutual 
diffusion 0,0]. 

This set of equations leads to two extended hydrodynamic diffusive modes. 

Since we are interested in exploring the packing aspects of jamming in a binary granular 
mixture which is controlled by the sizes rather than the masses of the particles, we confine 
ourselves to the case of equal masses mi = m 2 for the two species 151 ] . Here we analyze two 
illuminating special cases to bring out the role of size difference and packing in the jamming 
process which triggers the transition to a glassy state: (i) a system composed of a few large 
spheres in a matrix of small spheres and (ii) a system composed of a few small spheres in a 
matrix of large spheres. 

In both these extreme limits (x 2 = n 2 /n « 1 and X\ = rii/n << 1) the cross terms in 
the expressions representing the eigenvalues for the extended diffusive modes are negligibly 
small and the mode structure reduces to: 



(3) 



and 



z + (k) ~ -k 2 D(k)[J— f^k) + J^f 2 (k)} (4) 

V "1 V n 2 

Thus, there are two relevant modes: Z-(k), which governs the relaxation of total mass 
density fluctuations and z + (k), which governs the relaxation of concentration fluctuations. 
Let us analyze these modes for case (i). In this case, since the mixture consists of a large 
number of small spheres, the static compressibility Xr{k) which is the main determinant of 
the dynamics of a dense liquid, is given by Xr(k) ~ Xt^h(^)) i- e - ^ * s dominated by the 
static structure factor of the small spheres. Thus, z_{k) ~ ~^ 2 px o Sulk)-y L (k) • Consequently 
there is a significant slowing down of the dynamics of density fluctuations of the background 
matrix at the location of the first peak of Sii (A;) [isj] . In this case the static structure of large 



(type 2) spheres is flat and is given by S 22 (k) — 1. Thus, z + (k) ~ —k 2 Q D °, 2 k ) — —k 2 D, 
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where D 02 is the dilute gas limit of the diffusion coefficient of large spheres |1 71] . In case (ii), 
the large (type 2) and small (type 1) spheres switch roles. The final picture that emerges 
is the following. Caging and slowing down of dynamics is more effective for x 2 << 1 firstly 
because the ratio D G2 /D m = 0\jo 2 < 1 [3, Q] and a few large spheres diffuse slower in 
the background of small spheres compared to a few small spheres in the background of 
large spheres. In addition, in both cases there is also a significant slowing down of structural 
relaxation due to softening of the mode of total density fluctuations z_ (k) around the location 
of the peak of the static structure factor of the majority particles. 

In other words, caging is more efficient in a mixture dominated by small spheres compared 
to one dominated by large spheres. This is the main prediction made in this Letter for future 
experiments on dense granular mixtures designed to probe the efficiency of caging in such 
systems. 

The main point that we emphasize in this Letter is that by exploiting the analogy be- 
tween glass forming atomic liquids and granular matter we can draw some definite testable 
conclusions about caging dynamics in grains. 

To summarize, we have, for the first time applied extended or generalized hydrodynamic 
techniques to dynamics of granular matter. In particular, we consider a dense binary mixture 
of hard spheres and analyze the modes of density fluctuations which dominate the slow 
dynamics on length scales of the order of the sizes of the granular particles constituting the 
system. Our analysis points to some differences in jamming behavior between a mixture 



dominated by small spheres and one dominated by large spheres. Effects similar to the ones 



predicted here have been observed in confocal microscopy studies in colloidal mixtures [11]. 
The predictions made here can be tested against experiments in vibrated dense granular 
mixtures. We expect our analysis to be valid for a granular mixture consisting of nearly 
elastic spheres of comparable masses j^. While the present analysis captures the onset of 
glassy behavior in granular mixtures, it would be worthwhile to do a mode-coupling study 
for such a system using the extended hydrodynamic modes that stem out of our analysis 
as an input to understand glass transition in granular mixtures. In future, one can do a 
more complete analysis where the effects of momentum and temperature fluctuations and 
inelastic collisions [2l| are taken into consideration and check if the high density particle- 
scale dynamics presented here survives the inclusion of these effects. 
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